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In a previous paper ji]] an Euler angle parametrization for SU(4) was given. Here we present the 
derivation of a generalized Euler angle parametrization for SU(N). The formula for the calculation 
of the Haar measure for SU(N) as well as its relation to Marinov's volume formula for SU(N) 
H [| will also be derived. As an example of this parametrization's usefulness, the density matrix 
parametrization and invariant volume element for a qubit/qutrit, three qubit, and two three-state 
systems, also known as two qutrit systems will also be given. 
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I. INTRODUCTION 

The importance of group theory in understanding quantum mechanical processes has grown over the past 100 years 
as more and more physicists recognized that 

The basic principles of quantum mechanics seem to require the postulation of a Lie algebra of observables 
and a representation of this algebra by skew-Hermitian operators on a Hilbert space jq] . 

Today one finds applications of group theory in numerous research areas; from high-energy particle theory to exper- 
imental nano-scale physics. Any physical system that exhibits symmetries, or that can be thought of as Hermann 
indicated above, will have a group associated with it. In this paper we most concerned with representations of, and 
applications on, SU (N). This group shows up in numerous areas of study, most notably in numerical calculations con- 
cerning entanglement and other quantum information parameters. In order to assist in these numerical calculations, 
we have produced a parametrization of SU(N) that should offer some computational benefit. 

We will begin this paper by deriving a general Euler angle parametrization for SU(N). Afterward, a general 
equation for the differential volume element, otherwise known as the Haar measure, for SU(N) will be derived. Then 
we will show that this parametrization yields the familiar invariant volume element, generated from the integration of 
the Haar measure, for SU(N) as derived by Marinov [|| |j. Finally, the parametrization of N by N density matrices 
with regards to the general Euler angle parametrization for SU(N) will be shown. As an illustration of the usefulness 
of the parametrization, representations of qubit/qutrit, three qubit and two qutrit states will be given. 

II. REVIEW: EULER ANGLE PARAMETRIZATION FROM SU(2) TO SU(4) 

For a U £ SU(2) the Euler angle representation can be found in any good textbook on quantum mechanics or Lie 
algebras (see for example || ||) 

JJ _ e iC3Qi e iCT2Q2gio"3Ct3 (]\ 

For a U £ SU (3) the Euler angle parametrization was initially given in |)J Q and later in JlT[ 

JJ _ giA 3 ai giA 2 a2 giAsaa giA 5 a 4 ^iX 3 a 5 giA 2 c«6 giA 3 a 7 giA 8 a 8 fa\ 

For a U £ SU (4) the Euler angle parametrization was initially given in [Q 

JJ _ e iA3Qi e iA2Q2giA3a3 e iA5a4 e jA3a5 e iAioa6gjA3Q7 e iA2as 
x e iA3Q9 e iA 5 aio e i'A3Qn giA 2 ai2 giA 3 ai3 ^iXgau e iAi 5 ai5 

We would like to extend this work for a U £ SU (N) . 

III. LIE ALGEBRA FOR SU(N) 

From we already know how to construct the Aj for arbitrary SU(N). 1 
1. For every i, j = 1, 2, 3, . . . , N; i < j, we define two N by N matrices 

[A {2} (i, j)]^ = -i(S ifl 6j v - SwSjn), 
which form N(N — 1) linearly independent matrices. 2 



1 Georgi |12| also gives a method for constructing the N 2 — 1 Aj matrices for SU(N). 

2 We follow the standard physics practice by using hermitian generators. 
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2. Construct a further N — 1 matrices according to 
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for n = 2,3, ... ,N. 

By following this convention, N 2 — 1 traceless matrices can be generated. These matrices form a basis for the 
corresponding vector space and thus a representation of the SU(N) generators f|. For example, for N — 2 we 
generate the well-known Euler a matrices (i < j, < 2): 



n=[A {1} (1,2)]^ = 6 2fl 6 lv + 5 2u 6 



ifi ■ 



r 

1 or 



a 2 = [X {2} (1, 2)]^ = -i(S lfl 5 2v - SuthJ, 
and 



— i 

1 



<r 3 =A 



2 2 -l - 
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-1 
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2 2 - 2 VO -1 



(6) 



(7) 



(8) 



IV. DERIVING THE EULER ANGLE PARAMETRIZATION FOR SU(N) 

By following the work of Biedenharn jl3) and Hermann (5) we can now generate a Cartan decomposition of SU(N) 
for N > 2. 3 First, we look at the N by N, hermitian, traceless, matrices formulated in the previous section. This 
set is linearly independent and is the lowest dimensional faithful representation of the SU (N) Lie algebra. From these 
matrices we can then calculate their commutation relations 

[A,-, Aj-J = 2'/,,;, A/. . 

f ijk = ^Tr[[\ i ,X j ]X k ], (9) 

and by observation of the corresponding structure constants fijk one can see the relationship in the algebra that can 
help generate the Cartan decomposition of SU(N) (shown for SU(3) in || and for SU(4) in Q). 

Knowledge of the structure constants allows us to define two subspaces of the SU(N) group manifold hereafter 
known as K and P. From these subspaces, there corresponds two subsets of the Lie algebra of SU(N), L(K) and 
L(P), such that for ki,k 2 £ L(K) and pi,p 2 € L(P), 

[hi, k 2 ] e L(K), 
[pi.pa] eL(K), 

[k l7 p 2 ]eL(P). (10) 

For SU(N), L{K) = {Ai, . . . , A ( jv_i)2_i, Ajvj,} and L{P) = {A (A r_i)2, . . . , Aat2_ 2 }. 4 Given that we can decompose 
the SU (N) algebra into a semi-direct sum fl^] 

L(SU(N)) = L{K)®L(P), (11) 



3 The reason for this will become apparent in the following discussion. 

4 For example, for 5(7(4) we have L(K) = {X 1 ,X 2 , . . . , Ag, A i5 } and L(P) = {A 9 , A10, A14} |Q. By definition, L(K) and L{P) do not 
have any elements in common. This means that for N = 2, the second Ai element that is generated by this construction as an element 
of L(P) must be discarded. Similarly, the undefined Ao element in L(K) must also be discarded. For N > 3 one does not generate any 
duplications. 



4 



we therefore have a decomposition of the group, 



U = K P. 



(12) 



From pj we know that L(K) contains the generators of the SU(N — 1) subalgebra of SU(N), thus K will be the 
U(N — 1) subgroup obtained by exponentiating this subalgebra, {Ai, . . . , A(jv-i) 2 -i}i combined with Ajv2_i and thus 
can be written as (see [jjj ^| for examples) 



K = [SU(N-l)} -e iA « 2 -^ 



(13) 



where [SU(N — 1)] represents the (N — l) 2 — 1 term Euler angle representation of the SU(N — 1) subgroup. 

Now, as for P, of the 2(N — 1) elements in L(P) we chose the A 2 analogue, Xx su(N) , for SU(N) and write any 
element of P as 



P = K' ■ e iXx sv(.N)^ ■ K" 
where K' and K" are copies of K and ^x su{N) is given by the N by N matrix 



(14) 



= A(jv-i) 



/ 






\ i 



-i\ 






/ 



(15) 



NxN 



Unfortunately, at this point in our derivation, we have over parameterized U by 2N 2 — 6N + 5 elements 

U = K ■ K' ■ e iA («~i)2+i^ . K". 



(16) 



But, if we recall that U is a product of operators in SU(N), we can "remove the redundancies," i.e. the first K' 
component as well as the N — 1 Cartan subalgebra elements of SU(N) in the original K component, to arrive at the 
following product (again, see [0, H for examples) 



U = [SU(N — 1)] • e iA iv2-i^ e - iA iv2-i^ e _iA (iv--i)2~i0i e - iA (jv~2)2-i02 , . , g-iAa^JV-l 

x e <V-i)2 +1 <A . [SU(N - 1)] • e lA « 2 "i 
= [SU(N - 1)] • e -i V-i)»-i* • • • e^-'eVm+i* . [SU{N - 1)] • e iX » 2 -^ 



(17) 



By insisting that our parametrization must truthfully reproduce known vector and tensor transformations under 
SU(N), we can remove the last "redundancy," the final N 2 — 5N + 5 elements in K, 5 and, after rewriting the 
parameters, get: 



2<k<N 



^= II m)-[SU(N-l)}-e 



. pi* N 2_ 1 Ot N 2_ 1 



A(k) = e 4A3a (2^-3) e 4A ('=-i) 2 +i Q2 (' t - 1 ). (18) 
This equation, effectively a recurrence relation for the Euler angle decomposition of SU(N), can be further rewritten 



5 For N = 3, there are no redundancies. In fact, the -1, (3 2 —5*3 — 5 = —1), that occurs here is a result of removing one too many 
Cartan subalgebra elements from the end of K in the previous step. For TV = 3, one must restore a Cartan subalgebra element, in this 
case e 1 ^ 3 " 4 , back into its original position in the K component. For N > 3 this situation does not occur. 
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into a more explicit, and therefore final, form 

' " i n A ( k ^( N )) 



n A(k,j((N-i)) 

2<k<N ' v 2<fc<JV-l 

x e iA 3Q JV 2_ (JV _l) . . . g iA (JV-l)2-l Q iv2-2 e iA iv2-l Q N2_l 



Af>m>2 N 2<fe< 

x e iA 3a JV 2_ (JV _ 1) _ , , e iA (N-l)2-l Q i V 2_2 e iA iv2-l Q N2-l j 
A(k,j(m)) = e jA 3 Q (2fe-3) + J ( m ) e ^( fc -l)2 + l"2( fc -l) +J ( m)j 

fo m = N, 

j(m) = l J2 2(m + l) m^N. 

I 0<;<W-m-l 



(19) 



In this form, the important Xi matrices for equation (£L9h are 



(fc_l)2 + l 



A 



(JV-1) 2 + 1 



Ajy2- 




k < N, 







NxN 



-(iV-1), 



(20) 



NxN 



Notice that even though we restricted ourselves to N > 2 for the Cartan decomposition, equation ( |l9| ) is valid for 
N > 2. For our purposes it is enough to note that this parametrization is special unitary by construction and can be 
shown to cover the group by modifying the ranges that follow and substituting them into the parametrization of the 
characters |lT|]. 6 

V. PROCEDURE FOR CALCULATING THE HAAR MEASURE AND GROUP VOLUME FOR SU(N) 

Taking the Euler angle parametrization given by equation (|l9| ) we now develop the differential volume element, also 
known as the Haar measure, for the group. We proceed by using the method originally given in p!q ] and developed 
for SU{3) in §, |l| and for SU(4) in Q; take a generic U £ SU(N) and find the matrix 

U- l -dU (21) 

of left invariant one-forms, then take the determinant of the matrix of their expansion coefficients. 7 



6 One may be somewhat distressed by our "removal of the redundancies" statement that precluded the development of equation ( fis| ) so 
we offer another, geometrical, argument for the form of the parametrization in Appendix 

7 A detailed explanation of this method can be found in . We should note that one can take the determinant of the matrix of the 
expansion coefficients of the TV 2 — f right invariant one-forms which also yields the Haar measure in question. This is due to the fact 
that a compact simply-connected real Lie group has a bi-invariant measure, unique up to a constant factor. Such a group is usually 
referred to as 'unimodular' [H. 
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To begin, we take the transpose of U to generate 



u =U T 



[SU(N-l)] T -( [] A{k) T ) 



N>k>2 

J^N 2 -i a ff 2 -l e iX (N-l) 2 -i a N^-2 . . . e iX I a N 2-(N-l) 



X 

2<m<N x m>fc>2 



IJ LI A (k,3(m)) T ) (22) 



where 



A(k, j(m)) T = e 4A ^-D 2 +i Q2 < fc " 1 )+^'")e a 3«(2A ; -3) +3 ( m ) j 

f0 m = N, 

j(m) = l 2(m + l) m^N. ( 23 ) 



_0<i<Af-m-l 

An observation of the components of our Lie algebra sub-set 



{\fc-l) 2 + l, A fe 2_ 1 }, 

2 < k < N, (24) 



shows that the transpose operation is equivalent to making the following substitutions 

^(fe-i) 2 +l — * ~\fc-i) 2 +i, 

>&-i - ^-i, (25) 



for 2 < k < N in equation ( p2| ) generating 

U = e iA Jv2-i Q jv2_ le jA (N _ 1) 2_ 1 a JV 2_ 2 , . , e ^3a JV 2_ (N -i) 

x II ( II A (kJ(m))'), 

2<m<N \n>fc>2 ' 

A(k,j(m))' = e~ iA ( fc - 1 > 2 + lQ2<fc - 1)+J ' (m) e a3Q < 2fc - 3 >+ J '< m >, 

fo m = N, 

j(m) = l j2 2{m + l) m^N. ( 26 ) 

{0<l<N-m-l 

With this form in hand, we then take the partial derivative of u with respect to each of the N 2 — 1 parameters. In 
general, the differentiation will have the form 

= E L C( ai )E- L u (27) 

where we have used a "shorthand" representation with 

C(ai) e i * {— A( k _i)3 + i, A fc2 _ 1 }, 

2 < k < N, (28) 

and 

E L — e C(« N 2_,)a ]v 2_ 1 , , , e C(a !+ i)o!j + i 

= e -c( Q i+i) Q i+i . . . e -c( Q jv2-i) tt N2_i^ (29) 
By using these equations and the Baker-Campbell-Hausdorff relation, 

e x Ye- x =Y+[X,Y] + ^[X,[X,Y}} + ..., (30) 
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we are able to consecutively solve equation ( p7| ) for I — {N 2 — 1, . . . , 1}, leading to a set of N by N matrices which 
can be expanded in terms of the TV 2 — 1 transposed elements of the SU (N) Lie algebra with expansion coefficients cij 
given by trigonometric functions of the group parameters af. 



du 



M x ee ^-u- L = E L C{ ai )E- L = c y A J 



(31) 



The Aj's can be generated by using the methods contained in section III. 

Now the coefficients cy's are the elements of the determinant in question. They are found by evaluating the following 
trace [||: 



-Tr[\^-Mi\. 



(32) 



The index I corresponds to the specific a parameter, the j corresponds to the specific transposed element of the 
algebra. Both the I and j indices run from TV 2 — 1 to 1. The determinant of this TV 2 — 1 by iV 2 — 1 matrix yields 
the differential volume element, also known as the Haar measure for the group, dVsu(N) that, when integrated over 
the correct values for the ranges of the parameters and multiplied by a derivable normalization constant, yields the 
volume for the group. 

The full TV 2 — 1 by TV 2 — 1 determinant Det[cy], l,j S {A^ 2 — 1, . . . , 1}, can be done, or one can notice that the 
matrix can be written as 



a 



SU(N) 



c N 2 -l,N 2 -2 C N 2_ 1N 2_ 3 
C N 2 -2.N 2 ~2 c N 2 -2.N 2 -3 



C N 2 -1S C N 2_ 1N 2_ 1 
CjV 2 -l,l C N 2 -2,N 2 -1 



C\N 2 -2 



Cl,JV 2 -3 



Cl,l 



Cl,N 2 



(33) 



which differs only by an overall sign from Detfcy] above, but yields a quasi-block form that generates 

CSU(N) 






R 


T 


Q 



(34) 



In this form, R corresponds to the (N— l) 2 by (N— l) 2 matrix whose determinant is equivalent to dVi 



SU(N-l) 



■da 



N 2 -li 



Q is a complicated 2 (TV — 1) by (N — l) 2 matrix, and O is a (N — l) 2 by 2(N — 1) matrix whose elements are all zero. 
Now the interchange of two columns of a N by N matrix yields a change in sign of the corresponding determinant. 
But by moving 2 (TV — 1) columns at once, the sign of the determinant does not change, and thus one may generate a 
new matrix which is now block diagonal 



SU(N) 



R 





Q 


T 



(35) 



Thus, the full determinant is just equal to the determinant of the diagonal blocks, one of which is already known. So 
only the determinant of the 2(N — 1) by 2(N — I) sub-matrix T, 



Det[T] = 



C2(N-l),N 2 -2 ■■■ C 2 (N-1),{N-1) 2 
C\,N 2 -2 ■■■ Ci j (at_i)2 

is needed. Therefore the Haar measure for SU(N) is nothing more than 

dV SU (N) = Det[cjj] 

= -Det[T] * Det[D]da N 2_ 1 ...da x 

= -Det[T] * dV S u{N-i)da N 2_ 1 da 2 (N-i) 



(36) 



. da\. 



(37) 



This is determined up to normalization (explained in detail in Appendix |b|) . Integration over the 2 — 1 parameter 
space yields the following group volume formula (for A^ > 2) 



V f 



SU(N) 



dV< 



SU(N) 



dV< 



SU(N) 



(38) 
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where 

■»iV-2 



f^AT — 2 * ^(JV — 1) * N 

= 2 N - 2 * 2 N - 3 ■ ■ • 2° * 1 * 2 • • • (N - 1) * N 

= 2 e N\ (39) 



and 

e = £ (JV — I) 



2<1<N 

= (N - 2) + (N - 3) + (N - 4) + • • • + (N - (N - 1)) + (N - N) 
= (N - 2) + (N - 3) + (N - 4) + • • • + 1 + 

= i(AT-2)(iV-l). (40) 



VI. EXAMPLE: SU(5) PARAMETRIZATION, HAAR MEASURE, AND GROUP VOLUME 

Let U € SU(5). Following the methodology developed in the previous section, the Eulcr angle parametrization for 
SU(5) can be seen to be given by 



^ n ( n A (k,j(m))) 

5>m>2 x 2<fc<m ' 

x e iA 3 a 2 i e iA 8 a22 e iAi5a23 e iA24a24 

A(k,j(m)) = e iA3a < 2fc - 3 )+J<'")e jA ('=-i) 2 +i Q2 < fc - 1 )+^"*», 



m = 5, 

j(m) = <j £ 2(m + m^5. ( 41 ) 



0<K5-m-l 



Expansion yields 

u=( n ^,i(5)))( n ^( 4 )))( n ^(^h n ^(2)) 



x 2<fc<5 7 x 2<fc<4 7 x 2<fc<3 7 x 2<fc<2 

x gjA3Q21 gjA8Q22gjAi5Q23gjA24a24 

=A(2,j(5))A(3,j(5))A(4,j(5))^ 

^ giA 3 a 2 l giA 8 Q22 giAi s Q23 giA24«24 (42) 

The j( TO ) values are 

J (5) =0, 

i(4)= ^ 2(m + /) = 2m = 8, 

0</<5-4-l 

i(3)= 2(m + Z) = ^ 2(m + Z) = 2m + 2(m + l) = 14, 

o<;<5-3-i o<;<i 

j(2)= 2(m + /) = ^ 2(m + = 2m + 2(m + 1) + 2(m + 2) = 18, (43) 

CKK5-2-1 0<i<2 



and the A(k,j(m)) components are 



AO ifcW 


„iA 3 Qi iA 2 a 2 

— c o • 




AC\ i(^W 


„iA3Ct3 i\^an 

— o c , 




A(A 


„«A3ti5 i\ioae 

— CO , 




^W; J\°)) 


JA 3 a7 JAi 7 a s 

— o t; 3 




AO -i(A)) 


giA3ai +s iA 2 ci2+s _ 


. „iA 3 a:9 iA 2 Qio 

- o C , 




g«A3CK3-)_s giA5CK4 + 8 _ 


g*A3ttii g«A5Qi 2 


^(4,J(4)) 


= e iA3a5+8 e i ' > ' loQ6 + 8 




A(2,j(3)) 


giA3Qi + i4giA 2 a2+l4 


_ giA 3 ai 5 giA 2 Qi 6 


A(S,j(3)) 


giA3a 3+ i4giA5Q44.i4 


_ giA 3 Qi7giA 5 Qi 8 


A(2,j(2)) 


giA3ai + i 8 giA 2 a2+l8 


_ giA 3 Qi9giA 2 Q 20 



Thus 

jj _ e iA 3 Q!igiA 2 Q; 2 gjA3a3giA5Q4giA3Q5giAioQ!6giA3Q7giAi7a S giA3a9giA 2 QiOgiA3QiigjA5Qi 2 giA3Qi3giAioQi4 
x giA3ai5giA 2 ai 6 giA3Qi7giA5ai S giA3Qi9giA 2 a 20 giA3a 2 igiA8a 22 giAi5Q 2 3giA 2 4a 2 4 

with a differential volume element of 

dV S u{5) =~Det[T}dV SU ( 4 )da24,da$ ...dai 

=Det[T] cos(ai2) 3 cos(ai4) cos(ais) sin(2aio) sin(ai2) 
x sin(ai4) 5 sin(2cvi6) sin(ais) 3 sin(2a2o)cta!24 . . . da±, 

and where Det[T] is an 8 by 8 matrix composed of the following elements 

en ~Tr[X r Mi], 
— ■ i ( du 



(44) 



(45) 



(46) 



(1 < I < 8, 16 < j < 23, u = U T ), 



which when calculated yields 

Det[T] = cos(«4) 3 cos(ae) 5 cos(as) sin(2a2) sin(a4) sin(a 6 ) sin(as) 7 . 

Explicit calculation of the invariant volume element for SU(5) can be done by using equations (|3 
material from Appendix [b| From this one generates 

V su(5) = 2 6 * 5! * V(SU(5)/Z 6 ) 
72 ' 

which is in agreement with the volume obtained by Marinov j^]. 

It should be noted that in calculating Vsu(5) one must use the following ranges of integration from 
and expressed in equation ( |38| ) as V': 

< a 2l _i < 7T, < a 2i < ~, (l<i< 10) 

7T 

< OLl\ < 7T, < a 2 2 < —y=, 

v3 

7T 7T 

0<c*23<^=, < a 24 < — = • 
V6 v 10 



(47) 
(48) 

@ and the 

(49) 

Appendix H 



(50) 
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Note that these ranges do not cover the group SU(5), but rather SU(5)/Z^. The covering ranges for SU(5), following 
the work in Appendix are as follows 

< ai, Us, ai 5 , a w < 7r, 
< a 2 ,a 4: ,a e ,a 8 ,a w ,ai2,ai4,ai6,ai8 7 a2o < 
< a 3l a 5 , a r , a n , a i3 , a 17 , a 2 i < 27r, 
< a 22 < V37r, 

< a 23 < 2w|tt, 

< a 24 < y^TT. (51) 

VII. GENERALIZED DIFFERENTIAL VOLUME ELEMENT (HAAR MEASURE) FOR SU(N) 

Let us now quickly review the differential volume element kernels for the first few SU(N) groups (here rewritten 
in the order of the parametrization of the group): 

dVsu(2) = sin(2a 2 )da 3 da 2 dai, (52) 

dVsu(3) = sin(2a 2 ) cos(azi) sin(a4) 3 x dVsu(2)dasda4 . . . da\ 

= sin(2a 2 ) cos(a4) sin(a4) 3 sin(2a6)rfa;8 ■ ■ ■ dotx, (53) 

dVguU) = sin(2a 2 ) cos(a4) 3 sin(a4) cos(ag) sin(ae) 5 x dVgjj^da^daa . . . da\ 

= sin(2a 2 ) cos(a4) 3 sin(a 4 ) cos(ae) sin(ag) 5 sin(2ag) cos(aio) sin(aio) 3 sin(2ai 2 )(iai5 . . . dai, (54) 

dVsu(5) — sin(2a 2 ) cos(«4) 3 sin(a4) cos(ae) 5 sin(ae) cos(as) sin(as) 7 X dVsu(4)da24das ■ ■ ■ do>i 
= sin(2a 2 ) cos(«4) 3 sin(a 4 ) cos(ag) 5 sin(ag) cos(a§) sin(a 8 ) 7 

x sin(2aio) cos(ai 2 ) 3 sin(ai 2 ) cos(ai4) sin(ai4) 5 sin(2aig) cos(aig) sin(a 18 ) 3 sin(2a 2 o)do; 2 4 . . . da.\. (55) 

A pattern is emerging with regard to the trigonometric components of the differential volume element. For example, 
when one looks at the parametrization of the group and matches the trigonometric function in the differential volume 
element kernel with its corresponding e jA2 " m component, one sees sin(2a m ) terms showing up. It is plain to see 
then that, in general, the differential volume kernel is made up of trigonometric functions that correspond to group 
elements that are of the form e lA < fc - 1 ) 2 + 1 " 2( ' c_1)+; ' (m) and where j(m) is given in equation (|l9|). Therefore, there should 
be a general expression for the differential volume element kernel for SU (N) that can be written down via inspection 
of the Euler angle parametrization. We shall now show that this is indeed true, and give the methodology for writing 
down the differential volume element kernel for SU(N). We will also show that this procedure, after integration, 
yields Marinov's volume formula for SU (N) . 

To begin, let us take a look at the differential volume elements and their corresponding parametrizations for SU (3) 
and SU(4). From the parametrization originally given in g, |Io| for SU(3) we see that the e lAsQ4 term contributes 
to the cos(a4) sin(a4) 3 term in equation J53]). Yet in the parametrization for SU(4), originally given in |Q, the 
first e lA5Q4 term yields not cos(a4) sin(a 4 ) but rather cos (04) 3 sin (04). It is the second installment of A5 in the 
parametrization for SU(A), which can be seen to occur because of equation (|l8|), which gives us the cos(a4) sin(a4) 3 
term. For completeness, we should note that the cos(ae) sin(ag) 5 term in the differential volume element for SU(4) 
comes from the e lAloQ6 term in the parametrization. 

When we now look at the differential volume element and corresponding parametrization for SU(5), we see the 
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following relationships 



iA 5 a 4 _ 
c — 


COS^Cl4j SlIl^CK4j, 


6 — 


=> COS^g)^ sin(o!g), 


^i\nas 


cos(as) sin(ag) , 


giA5Q:i2 


cos(ai 2 ) 3 sin(ai 2 ), 


giAioaw 


=4> cos(ai4) sin(ai4) 5 , 


giAsais 


cos(ai8) sin(ais) 3 . 



(56) 

By combining all these observations from the parametrizations and differential volume element kernels for SU(2) to 
SU(5) we can see that the following pattern is evident 



o iA (N-l)2 + l Q 2(JV-l) 



sin(2a 2 ) 

cos(a 2 ( fc -i)) 2fc ~ 3 sin(a 2 (fe-i)) 
cos(a 2 (jv-i)) sin(a 2(A r_ 1) ) 2Ar ~ 3 



k = 2, 
2 < k < N, 
k = N, 



(57) 



relating the A(k) term of the recurrence relation given in equation ( |18[ ) with the Det[T] term of the differential volume 
element given in equation (p7|) . Exploiting this recurrence relation we now have a methodology for writing down the 
kernel for any SU(N) differential volume element. Recalling equation (jn|) and the A(k,j(m)) term we see that in 
general 



e iA (fc-l)2+l Q 2(fc-l)+j( m ) 



sin(2a 2 +j(m)) 

c os(a 2 (fc_i) +:) ( m )) 2fc_3 sin(a 2 (fc_i)4.j( m )) 
cos(a 2 ( m _i)+j( m )) 8m(a2( m -i)+j(m)) m ~ 



k = 2, 
2 < k < m, 
k = m, 



(58) 



for N > m > 2. Using this result, combined with the knowledge that only these parameters contribute to the 
integrated kernel in dVsu(N) , we are able to write the following product relation for the kernel of the Haar measure 
for SU(N) 



where 



K 



SU(N) 



Ker(k, j(m)) 
and j(m) is from equation (jlo|) . 



JJ [] Ker(k,j(m))), 

N>m>2 ^2<k<m 

sin(2a 2+j(m )) 



. dai 



(59) 



cos(a 2 ( fe _ 1 ) +J („ l )) 2fe 3 sin(Q!2(fc_i) + j( m )j 
cos(a 2 (m-i)+j(m)) sin(a 2 ( m _i) +J ( m )) 2m_3 



k = 2, 

2 < k < to, 

fc = TO, 



(60) 



VIII. EXAMPLE: SU(6) HAAR MEASURE AND VOLUME CALCULATION 

As proof of the validity of equation (|6(]) we shall use it to write down the differential volume element for SU(6). 
Observation of equation (|37]) tells us that dVsu(6) 1S dependent on the differential volume elements of SU(5), St/ (4), 
SU(3), and SU(2). Thus, in the process of writing down dVsu{&) we w iU n °t only confirm the calculated differential 
volume elements for the previous four SU (N) (TV = 5, 4, 3, 2) groups, but we will also be able to give the functional 
form of Dct[T] in equation ( |37| ) for 5/7(6) without having to formally calculate the 10 by 10 determinant. So, for 
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N = 6, equations (Q) and (|T^) yield 

6>m>2 ^2<fc<m 



sin(2a 2+:) („ l )) k = 2, 

Ker(kJ(m)) = { cos(a 2{k _ 1)+j{m) ) 2k ' 3 sm(a 2(k _ 1)+:j{m) ) 2 < k < 

(«2(m-l)+j(m)) s i n («2(m-l)+j(m)) 2m ~ 3 k = m, 



m. 



cos 



m = 6, 

j{m)=( £ 2(m + J) m^6. ( 61 ) 



0<i<6-m-l 



Which when expand gives 

^St/(6) = ( II ^er(fc,j(6))) ( [] #er(fc,j(5))) ( JJ tfer(fc,i(4)) 

2<fc<6 ' ^2<fc<5 ' 

JJ tfer(fc,j(3))) C [] Ker(k, j(2))) 



x 2</c<3 7 v 2<fc<2 7 

= (Ker(2, j(6))#er(3, j(6))Ker(4, j(6))Ker(5, j(6))Ker(6, j(6))) 
x (Ker(2, j(5))Ker{3, j(5))Aer(4, j(5))Aer(5, i(5))) 
x (Xer(2, j (4))if er(3, j(4))Aer(4, j(4))) 
x (Ker(2,j(3))Ker(3,j(3))) 

x (Ker(2,j(2)). (62) 



The j(m) values are 
j(6) -0, 

i(5)= XI 2(m + /) = 2m = 10, 



0<K6-5-l 



j(4)= X 2(m + /)= X 2(m + i) =2m + 2(m+l) = 18, 

0<i<6-4-l 0<i<l 

j(3) = X 2(m + Z)= V] 2(m + J) = 2m + 2(m + 1) + 2(m + 2) =24, 

0<;<6-3-l 0<i<2 

j(2)= 2(m + /)= 53 2 ( m + =2m + 2(m+l)+2(m + 2) + 2(m + 3) = 28, (63) 



0<i<6-2-l 0<2<3 
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and the Ker(k, j(m)) components are 






Ker(2,j(Q) 


= sin(2a 2 ), 




Ker(3, 7(6) 


= cos(a4) 3 sin(a4), 




Kev(4: 

V 7 ./ V / 


= cosf ttfi)^ sin(Qfi) 




Kev(5 i(o) 


= cos(cks)^ sin(o;8), 




Kev(G i(o) 


v< V *- ± u / " *v J-U / 7 




Aer(2, 7(5) 


= sin(2a?^m) = sin(2a;-|9), 




Ker(3 7(5) 


= COSfc^^inl sinfcK/i-Li nl = 


cos(o;i4) 3 sin(o;i4) 




= COsfo^-Li n 1 Sill fr>K_Li nl = 


oosfo'-ifi r sinfn'if;') 




— cnsT rvo 1 ml si n frvo i 1 nl^ — 


pnsfrvi o 1 sin Trv-i o 1 ^ 


Ker(2 


= sin ( 2o , oxi « ) = sin(2o , on) 




Ker(3 i(4) 




msf'rvon 1 stin (rvon 1 
lajo ^Lt 2 J } oaii^uey^ y 


7>" /I ' / /I \ 

A er(4, j (4 J 


= cos(a 6+ i 8 ) sin(a 6 +i8 ) = 


„„„/ \ • / \5 

cos(a 2 4j sm(a 2 4j 


Ker(2,j(3) 


= sin(2a 2+24 ) = sin(2a 26 ), 




i\Ter(3, j(3) 


= cos(a 4+24 ) sin(<xi +24 ) 3 = 


cos(a 28 )sin(a; 28 ) 3 


*fer(2, j(2) 


= sin(2a 2+28 ) = sin(2a 30 ). 





Thus 



A 



S(7(6) 



: sin(2a 2 ) cos(o!4 

x sin(2ai 2 ) cos 

x sin(2a!2o) cos 

x sin(2a26) cos 
x sin(2a 30 ) 



(64) 



3 sin(a4) cos(«6) 5 sin(ae) cos(as) 7 sin(ag) cos(aio) sin(aio) 9 
an) 3 sin(«i4) cos(aie) 5 sin(aig) cos(ai8) sin(aig) 7 
a 22 ) 3 sin(a 2 2) cos(a 2 4) sin(a 2 4) 5 
a 28 )sin(a 28 ) 3 



and 



SU{6) 



Ksu(6)da 3 



. da±. 



(65) 
(66) 



Comparison of the above kernel with those from equations (|52|), (|53"|), (|54|), and ( p5| ) confirms that equation ( |60| ) does 
correctly yield the differential volume element for SU(6). As added proof, integration of Ksu(6) using equations (|3^ ) 
and (J39|) combined with the following ranges (the general derivation of which can be found in Appendix 0) 



< a 2j _i < 7r, 
< a 3 i < 7T, 



< a 33 < 



< a 35 < 



V6' 

7T 



< a 24 < -, 
0<o 32 <^, 

< a 34 < 



(1 < i < 15) 
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15 



(67) 



yields 



V f 



SU(6) 



,10 
7T" 



2 lu *6!*y(S'f/(6)/Z 6 ) 

r 20 



1440V3 

which is in agreement with the volume obtained by Marinov 



(68) 



IX. GENERALIZED GROUP VOLUME FOR SU(N) 



In looking at equations (|38j) and ( J60| ) we see that under this Euler angle parametrization for SU(N) the calculation 
of the invariant group volume is simply a matter of successive integrations of sines and cosines, multiplied by some 
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power of tt and a normalization constant. Therefore it stands to reason that with the derivation of a generalized form 
for the differential volume element of SU(N), i. e. the Haar measure, for SU(N), there should be a corresponding 
generalized form for the volume element for SU (N) . It is to this derivation that we now focus our attention. 

We begin by noticing that in our Eulcr angle parametrization we have a total of N 2 — 1 parameters, of which the 
final N — 1 are the Cartan subalgebra elements for SU(N), thus leaving N(N — 1) elements evenly split between the 
-\fc-i) 2 +i an d A3 parameters (for 2 < k < N). Rewriting this observation using equation (19) we see 

N(N - 1) »V_ 1)2+1 a 2 (fc-l)+3Cm) 



2 

N(N — 1) 



e 



2 

TV - 1 => e lA 3" Ar2 -( Ar - 1 ) _ _ , e iA N2-i«jv2-i^ (69) 

Examination of equations ( |59| ) and ( |60| ) shows that the ./V — 1 Cartan subalgebra elements and the N(N — l)/2 A3 
elements do not contribute to the integrated kernel, but their corresponding parameters are integrated over. Expanding 
on the general results from Appendix [b|, we see that we can use the following ranges to calculate the group volume 
for SU{N) 

tt N(N-l) 
< a 2i _i < tt, < a 2l < -, 1 < 1 < '-, (70) 

and 



0< a N 2 +b <ttJ—, —r , a ee b mod (N + 1) a = 2...N. (71) 

V a ( a - !) 

From these ranges, it becomes apparent that each A3 element contributes a factor of tt to the total integration of the 
differential volume element over the iV 2 — 1 parameter space while each of the N — 1 Cartan subalgebra elements 
contributes not only a factor of tt but a multiplicative constant as well. Explicitly 

N(N-1) 

da 2 i-i — tt 1 < 1 < 

W(jV-l) 

==>• tt 2 from the A3 components, (72) 



and 



Jo 



da^[2 +b — tt * — a = 6 mod (TV + 1) a = 2...N 



a(a — 1) 



fc(fc-l) 



2 < k < N 



IT % i — from the N-l Cartan subalgebra components. (73) 

2 <i<iv V^- 1 ) 

We now focus our attention on the integration of the differential volume element kernel given in equation (|6C|). 
Examination of the Ker(k, j(m)) term combined with the previously given ranges, yields the following three integrals 
to be evaluated 

sm(2a 2+ j( m ))da 2+ j( m ) k = 2, 

coa(a 2 (k-i)+j( m )) 2k ^ 3 sm(a 2 (k-i)+j{m)))da 2 (k-i)+j( m ) 2 < fc < m, 

cos(a 2 ( m -i) +:) („ l ))sin(a2( m _i) + j( m )) 2 ' n ~ 3 )(iQ!2(m-i)+j(m) k = m, (74) 
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and where, again, j(m) is from equation (|l9|). The first integral is equal to 1, and the other two can be solved by 
using the following integral solution from Dwight's "Tables of Integrals and Other Mathematical Data" 



sin p (x) cos q (x)dx = nT J pH _ ? - \ / p+l,g + l>0 (75) 



r(£±i)r(a±i 
2r(4^ + i) 

where T(x) is the standard Gamma function with the following properties 



r(i) = i, 

r(ar+l) =a:r(x). (76) 



Thus when p — 1 and q — 2k — 3 we get 



r( i±i )r( 2fc^+i ) ^ r(i)r(fe-i) 
2 r(i+f^ + 1) ~ 2r(fc) 

r(fe-i) 



2(fc - i)r(fe - 1) 
i 



2(fc- 1) 



(77) 



and when p = 2m — 3 and g = 1 we get 

r( 2m -3+i ) r (i±i) ^ r(m- i)r(i) 

2r( 2m - 3+1 + 1) ~ 2T(m) 

r(m-l) 



2(m- l)r(m-l) 
1 



2(m- 1)' 



(78) 



From this work we can see that with regards to the integration of the Ker(k, j(m)) terms in equation ( |60| ) their 
contribution to the overall calculated volume clement is 

/ II ( II Ker(k,j(m)))da 2i = U ( U V ( fc < m ) 

J ° N>m>2 ^2<fc<m ' N>m>2 ^2<fc<m 



and where 



NlN - 1) 

dan runs from 1 < i < , (79) 



1 k = 2, 

nKm) = \ , 2<k<m _ (80) 



I 2(fc-l) 

Substituting this result as well as equations J7^), (|73]), and (|5|) into equation (H|) yields 
^S(7(Af) = / ■ • • y ^sc/(Af) 



Q N * J ■■■ J K SU ( N )da N 2_ 1 ...dai 
v 



(JV-l)(jV-2) N(N-l) 



n ^m^t) n n *»,».>)). w 

2<k<N V v y/ x AT>m>2 v 2<fe<m 
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This expression can be simplified by observing that 



n 'feV 2 ^'' n Jj^y, 

2<k<N V v ' 2<k<N V v ; 



1 1 1 

* ; * ; * • * ■ * 



V2 * 1 V3*2 V4*3 a/TV * (AT - 1) 



2^^-x) 1 ( 8 2) 

y/N\(N- 1)! 



and, through the usage of equation (p0|), that 



n ( n v(*,m)) = ( n v^^f n m^-i)) 

>m>2 ^2<fc<m ^2<k<N ' ^2<k<N-l ' 

n v (M))( n v(fe,3))v(2,2) 



X • • • X 

^2<fc<4 7 v 2<fc<3 

=(V(2, AT)V(3, AT) • ■ ■ V(AT, A r ))(V(2, N - 1)V(3, AT - 1) • • • V(AT - 1, N - 1)) 

x ■•• x (V(2,4)V(3,4)V(4,4))(V(2,3)V(3,3))V(2,2)) 
=V(3, AT)V(3, AT - 1) • • ■ V(3, 4)V(3, 3) x V(4, AT)V(4, N - 1) • • • V(4, 4) 
x V(AT - 1, N)V(N - 1, AT - 1) X V(AT, N) 
l\ N ~ 2 fl\ N ~ 3 ( 1 \ 2 / 1 



X - X • • • X — — — X 



4/ \6/ \2(N~2)J K 2 ^- 1 ) 

.\N-2,^\N-2 fl\ N ~ 3 /l\ N ~ 3 /l\ 2 / 1 \ 2 /l\ / 1 



X - - X •• • X 



2) \2 J \2J \3J \2 J \(N-2)J \2 J \(N - 1) 



N-l 

(N—2)(N—1) 

2 s 



fe=i v ' 

By substituting these results back into equation (|8l"|) we get 



N-l 



W-l (JV-l)(iV + 2) 



y/N\(N- 1)1 



k=l x 

/ 1 \ 

which is just Marinov's initial formulation of the volume of SU(N) [|. This is an important result, for it shows 
that the overall generalized Euler angle parametrization of SU (N) gives results that are consistent with previously 
scrutinized work which used a completely different methodology to derive the invariant volume of SU(N). 8 

X. N BY N DENSITY MATRIX PARAMETRIZATION 

We now turn out attention to the parametrization of N by N density matrices. We state that one may parameterize 
any N hy N density matrix as (see |T3] for more details) 



p = U Pd U\ (85) 



Marinov calculated the invariant volumes for all compact simple Lie groups by exploiting the spectral expansion of the Green's function 
for diffusion on a group manifold M. 
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where pd is the diagonalized density matrix which corresponds to the eigenvalues of the (N — l)-sphere, S 



N-l 



cos 2 (6>i)---sin 2 (6>jv-i) 



V 











^ 



COS 2 (6»7V-l)/ 



NxN 



Mi, 


o . 


. 





A 2 . 


. 


V o 


. 


■ A TV 



(86) 



NxN 



the range of the 9 parameters is given by J!?], [T^] 



cos 



1 



< h < -■ 
2 



(87) 



and U is from equation ( jl9|). For completeness we should note that W is defined to be the conjugate transpose of U, 
which, through equation (pq), can be seen to be equal to 

Jji = e ~ iA lv2-l Q N2_l e - iA (JV-l)2_l Q iv2-2 . . . g _i A 3 a N? - (N - 1) 



n ( n Mk,j(m)y), 



2<m<N x m>fc>2 

A(k, j(m)) f = e _4A ( fc - 1 ) 2 + lQ2<fc - 1 '+ j(m) e- a3Q < 2fc - 3 '+^-», 

{0 m = N, 

J2 2(m + m^iV. 
0<Z<A r -m-l 



Throughout the rest of the paper, will be parameterized by the following set of quantities |l], [Ti| [IT], 18 



2<n<N 



f(9i,02,---,0N-i) = -Tr[p d ■ Xn 2 -l\- 
For example, for the density matrix of two qubits U is given by equation (|3|), is 

jj\ _ p — iAi5Qi5 g— iXgau^— 1X30113 g— iA2«i2 g— iAsaii g— iAsaiOg— iAsag g—iA2 cte 



x g— iA3Q7g — iAioa6g— iA3Q5g— iA 5 a 4 g— iA 3 Q 3 g — iA 2 Q2g— iA 3 ai 



and is given by 



Pd 



/sin 2 (0i)sin 2 (0 2 )sin 2 ( 



V 



\ 

cos 2 (e»i)sin 2 (6»2)sin 2 (6l 3 ) 

cos 2 (6> 2 ) sin 2 (#3) 



'3) 

at a \ „:„2^ a \ „ ; „2 











where 



\U + + 2uj 2 )x 2 y 2 * A 3 + -j= (-2 + 3x 2 )y 2 * A 8 + ^=(-3 + 4y 2 ) * A 15 



w 2 = sin 2 (ex), x 2 = sin 2 (0 2 ) , y 2 = sin 2 (0 3 ) , 



7T „ 7T 
- < 01 < -. 



COS _1 (-F= )< #2 < -, -<0 3 <-, 

V V3 ~ 2' 3 _ 2 



(89) 



(90) 



(91) 



(92) 



and the one-quarter normalization of I4 keeps the trace of pd in this form still unity jl). Therefore, using equations 
(|l9|), (|35|), (88), and (|B9"|), one can easily write down any N by AT density matrix. 



18 



XI. EXAMPLE: HAAR MEASURES, GROUP VOLUMES AND DENSITY MATRICES FOR 
QUBIT /QUTRIT, THREE QUBIT AND TWO QUTRIT STATES 

Using the formalism we have now established, it is quite easy to write down the Haar measure and group volume 
for ST/ (8) and SU(9), as well as the 6 by 6, 8 by 8 and 9 by 9 dimensional density matrices, that correspond to the 
qubit/qutrit, three qubit and two qutrit states. 9 

Case 1, Qubit/ 'Qutrit states 

Here, a qubit interacts with a qutrit (a three-state system), thus yielding a 6-dimensional Hilbert space which needs 
to be parameterized. Using equations (|l9|), ([Bo]), (|8q), and (|89|) we arrive at the following formula for the density 
matrix of a qubit/qutrit system: 

u = n ( n A (k,j(m))) 

6>m>2 ^2<fc<m ' 

x giA 3 a 3 i gjA 8 a32 giAi 5 a 33 g2A 24 a 34 giA 35 a 3 5 

A(k,j(m)) = e lA3Q ( 2fc - 3 )+^™)e iA ( fc - 1 ) 2 + ia2 < fc - 1 »+ J < m », 

Jj\ _ g— iA35CK35 g— iA24«34g— iAl5a33g— i^8a32g— iA3«31 



Explicitly: 



n ( n Mkj(m)y), 



X 

2<m<6 \n>k>2 



{0 m = 6, 

Y, 2(m + l) m^6, 
0<i<6-ro-l 

Pd=^6+ A)*A n2 _!. (93) 

2<n<6 



„ _ giA 3 Ql gjA2Q2 giA3Q3 giA5Q4giA 3 Q!5 giAl0Q6 giA 3 a7p«Al7as 

x giA3Q9 e iA26aiOgiA3aii e iA2Qi2giA3ai3 e iA5ai4 e iA 3 ai5 e iAioQ;i6 

x giA3ai7 e iAi7ai Se iA3ai9 e iA 2 a 2 OgiA3Q!2lgiA5a22giA3Q!23gi- > 'loa 2 4 

x giA 3 a 2 5gi- > ' 2 a26giA3a27giA5a2SgiA3a29giA2Q30 (p^ e ~ »^203Q g — iA 3 a 29 

v _ — iA5«2S „ — iA3Q27 „ — iA2Q26 „ — iA 3 C(25 „ — iAloa24 „ — iA 3 C(23 „ — iAsQ22 _ — iA3C(21 

x g-iA2a20g— *A 3 ai9g— iAi7ai 8 g— iA 3 ai7g— iA 10 ai6g— iA 3 ai5g— iAsaiig— iA 3 ai3 
x g — iA 2 ai2g — »A 3 aii g — iA 2 6aaOg— iA 3 a 9e — iA^agg— iA 3 Q7 e — iAioQ6g — iA 3 a 5 

x e^ lA:sQ4 e~ lA3a3 e- iA2Q2 e- iA3Ql , (94) 



The Haar measure and group volume for SU (6) has already been written down and calculated in section VII] 
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where 



1 cos(2 6li)sin(6l 2 ) 2 sin(6'3) 2 sin(6'4) 2 sin(6'5) 2 
Pd = g 1 6 2 * 3 

(1 + 3 cos(2 6 2 )) sin(6'3) 2 siri(6> 4 ) 2 sin(6> 5 ) 2 

* As 



4\/3 

(1 + 2 cos(2 6> 3 )) sin(6> 4 ) z sin(6» 5 ) 

2~7i 

(3 + 5 cos(2 4 )) sin(0 5 ) 2 



4V10 
2 + 3 cos(2 6> 5 ) 



* A 



24 



(95) 

From Appendix [§] and equation (^7|) we know that the ranges for the quotient group SU(6)/Zq and the pd parameters 
are 

< ctM-i < vr, < a 2i < |, (1 < i < 15) 

cos" 1 f ^=y) < 9j < \ (1 < 3 < 5)- (96) 

And, for completeness, we again note that the Haar measure and group volume for a qubit/qutrit system has already 
been calculated in equation (|68|). 

Case 2, Three qubit states 

Here, three qubits interact, thus yielding a 8-dimcnsional Hilbcrt space which needs to be parameterized. Using 
equations @, @, @, and @ we arrive at the following formula for the density matrix of the three qubit system: 



i7 = n ( n A (kj(m))) 



8>m>2 x 2</c<m 
x giA 3 a57 g»AsQ:58 giAi 5 Q59 giA 2 4a60 giA3 5 a 6 l giA 48 a 6 2 giA 6 3C«63 

A(fc,j(m)) = e lA 3«(a-3) +J W e '\ )i -l)2 + l«2( t -l) +I ( m ) j 



n ( n A{k,j{m)A, 



X 

2<m<8 v m>fc>2 



A(kJ(m)^ = e -^(fc-i)= + i a 2C=-i)+i(™) e -iA3a (afc _ 3)+J(m , s 

[O m = 8, 



0<2<8-m-l 



4 ls+ 51 f(0i,0 2 ,...,e 7 )*X n ,_ 1 . (97) 



2<n<8 
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Explicitly: 



giA3a7gjAi7Q8 e iA 3 Q9 e iA26ai0g^3ail gjA 37 ai2 



X e 

A2Q16giA3Q17giA5ai8 



^ g«A3Q13 giA5oai4 giA3Q15g^v ZLllD g t 
x gjA3Qi9glAi Q20g«- > »3a21 g « A17 a 2 2 g « A3 Q 2 3 g^26 "24 

x giA 3 Q 2 5 gjA3 7 Q26giA3a27gi^2Q;28 giA3C«29 gj->»5a30 

w c iA3a3l iAl0Q32 o iA3Q33 ^iAl7a34 2A3a35 „?A26«36 
AC c- C c- C c- 

x g^3"37 g«A 2 Q3 8 gjA 3 Q3 9 giA 5 Q 40 giA 3 Q 4 i giA 10 a 42 

x e i^3°"t3 gi^naa giA 3 a45 g* A2«46 gjA 3 a 4 7 giA 5 a 48 

x ^^3^49 g^lO^O giA 3 Q5i glA 2 Q52 giA3Q 53 giA 5 Q 5 4 

x g»>'3a55gj>'2a56 (^p d ^ e -i^20i56 g-«A3«55 g-iA 5 a 5 4 g-iA 3 Q 5 3 

x g — iA 2 a52g — iA 3 a 5 i g — iAioasOg — 'iA3a 49 g — iA 5 a48g — '^3047 

x g — iA 2 Q46g — jA 3 Q45g — iAi7Q 44 g — iA 3 Q43g — iAioQ 4 2g — iA 3 a 4 i 

x g — « A 5 Q 40 g — iA 3 a39 g — iA 2 Q38 g — iA 3 a 3 7 g — iA 2 6a36 g — ^30:35 

x g — «Al7Q34g — jA 3 a33g — iAioQ 3 2 g — iA 3 Q 3 i g — jA 5 Q30g — lA 3 a29 
x g — «A2Q28g — jA 3 a27g — «A37Q2Bg — iA 3 a 2 5 g — iA26CK24g — jA 3 a 2 3 
x g — »Ai7Q 2 2 g— iA 3 Q21 g — iAioQ 2 Og— iA3Qi9g — lA 5 Qi8giA3Qi7 

4 g— t«a«ia g — *A37Qi2 g — 'iA3Qn 



x e 
x e 



— lA 2 Ql6g — iA 3 Qi5 g — iA 5 oQi4g — iA 3 Qi3g 



x e 



— iA 2 6Ql0g — 'iA3Q 9 g — iAi7Q 8 g— iA3Q7g— iAi a6g— iA 3 a 5 

— iA 5 Q 4 — iA 3 Q3 — iA 2 a2g— iA 3 Qi 



(98) 



where 



Pd = 



1 8 cos(2 6li)sin(6'2) 2 sm(6'3) 2 sin(6'4) 2 sin(6' 5 ) 2 sin(6i 6 ) 2 sin(6'7) 2 . 

*A 

8 2 

(1 + 3 cos(2 9 2 )) sin(<9 3 ) 2 sin(<9 4 ) 2 sin(<9 5 ) 2 sin(0 6 ) 2 sin(6» 7 ) 2 x 

7= 

4^3 



(1 + 2 cos(2 6> 3 ) 



(3 + 5 cos(2 6> 4 ) 



_ (2 + 3 cos(2 6> 5 ) 



(5 + 7cos(2 6> 6 ) 



4V2T 

(3 + 4 cos(2 6 7 ) 

4~77 



4^3 

3in(6*4) 2 sin(6' 5 ) 2 sin(# 6 ) 2 sin(# 7 ) 2 
= *A 

9 



2^ 

sin(6^ 2 sin(6» 6 ) 2 sin(6» 7 ) 2 

7= *A24 

4^ 

sm(6> 6 ) 2 sin(6> 7 ) 2 

7= *^35 



■15 



T5 
sm(d 7 ) 



*A48 



*A 6 3- 



(99) 



From Appendix [b] and equation (^7|) we know that the ranges for the quotient group SU(8)/Zg, and the pd parameters 
are 



< a2i-l < 7T, 



< a 2i < |, (1 < i < 28) 



cos 



-1 



(1 < J < 7), 



(100) 
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combined with equations ([59]) and (|6C)|), which yields the Haar measure for SU(8) 

dVsu(s) = sin(2a2) cos(a4) 3 sin(o!4) cos(ae) 5 sin(ae) cos(as) 7 sin(ag) cos(aio) 9 
x sin(aio) cos(ai2) 11 sin(ai2) cos(ai4) sin(ai4) 13 

x sin(2ai 6 ) cos(a i8 ) 3 sin(ai 8 ) cos(a 2 o) 5 sin(a2o) cos(a 22 ) 7 sin(a 2 2) cos(a 2 4) 9 
x sin(a 2 4) cos(a 2 6) sin(a26) 11 

x sin(2a 2 8) cos(a 30 ) 3 sin(a 3 o) cos(a 32 ) 5 sin(a 32 ) cos(a 34 ) 7 sin(a 34 ) cos(a 36 ) sin(a 36 ) 9 
x sin(2a 3 8) cos(a4o) 3 sin(a4o) cos(tt42) 5 sin(a42) cos(«44) sin(a44) 7 
x sin(2a46) cos(a4s) 3 sin(a4g) cos(a5o) sin(aso) 5 
x sin(2a52) cos(as4) sin(a54) 3 

x sin(2a56)(ia6 3 . . . dai, (101) 
one can calculate the group volume for the three qubit system; 



8-1 (8-l)(8+2) f-f-r I 1 \ 

^(8) = 2 ^^^^ II (fell 

fc=l ^ ' J 



391910400 

which is what one gets with N — 8 from equation (p3). 



(102) 



Case 3, Two Qutrits 

Here, two qutrits (two three-state systems) interact, thus yielding a 9-dimensional Hilbert space which needs to be 
parameterized. Using equations (^), (^5|), (|S8|), and ( |89| ) we arrive at the following formula for the density matrix 
for two qutrits: 

U = 



[] ( II A (k,j(m))) 



9>m>2 x 2<fc<m 

w A\zOLf& zAsa74,p*Al5CK75 i\lA,Otl& Jl\zh<^ll „iA48«78 „i)"QZOL79 AXsqolsq 
A c fc. o fc. o fc- fc- fc- , 

A(k,j(m)) = e jA3Q ( 2 ^ 3 )+^™)e lA < fc - 1 ) 2 + ia2 < fc - 1 >+ J <"*», 

Jj\ _ g — iAgQQsO g — 4^63 «79 g~ iA 4s a7 S g — i A35 a 77 g — iA 2 4Ct76 g~ 4^150:75 g — i\ s Q74 g— iAa 0:73 



n f n A(A,i(m))tV 



x 

2<m<9 v m>fe>2 



A(fc,j(m)) 1 ' = e^^C 1 - 1 )^ 1 " 2 ^ -1 '^'" 1 ' e~ iA3a ( 2fc - 3 )+3( m ) , 

fo m = 9, 

j(m) = .j 2(m + m^9, 

I 0<Z<9-m-l 



~1 9 + ^ f(9u6 2 ,...,0 S )*K*-i- (103) 



2<n<9 
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Explicitly: 



where 



p _ giA3ai giA2«2 giA3a3 giA5Q4 giAsas 



x giA 3 a7 e iAi7as giA 3 a 9 giA26aiOgiA 3 Q!ii g iA 



giA26ai0giA 3 Q!iigiA37ai2 
x giA 3 ai 3 giA5oai4giA 3 Q!i5 ( 



5giA 6 5ai 6 giA 3 ai7giA2Q!is 
x giA3ai9giA5C«20giA3Q2lgiAioo 



1 giAl0Cl22 giA 3 C«2 3 giAl7024 

x giA 3 a 2 5giA 2 



gi^ 3 a27g*A370 



x g»^ 3 "25giA 2 6Q;26gj- > i3a27giA37Q28giA3Q:29giA5oa30 
x giA 3 Q!3igiA2Q!32giA3Q33 £ 
iA-jrv-j-r iAi-7fv^o jA^ftsn 



r g<>^37 <-*2K gt^3Lt2y g'^5U<^; 
3 giA 5 Q!34giA3Q!35 giAi Q!36 

x giA 3 Q!37giAi7a3 8 giA 3 0!39 giA 2 eQ;40 giA 3 0!4i giA 3 7a 4 2 

x giA 3 Q!43 giA 2 0!44giA3Q!45 giA 5 a 46 giA 3 0!47giAioa!48 

x giA 3 Q!49 giAi7a so giA 3 Q!5i giA 2 ea52 giA 3 0!53 giA 2 054 
x giA3a55giA 5 Q!5 6 giA3a57giAi 
x giA3a 6 l „i>> 



giAioass giA3a 59 giAi7a 60 

giA 2 Q!62giA3 



x giA 3 0!|37 iA2C«68giA 



: giA3Q! 6 3 giA 5 a 64 giA3Q!i35 giAioaee 
iA 3 Q! 6 9 iAsi 



giA 5 O!70giA3Q!7igiA2Q!72 

>2Q;72g— iA 3 a7i g — iA 5 O!70g— iA 3 0! 6 9 

■3"67g — iAi Q!66 g — iA 3 0!|35 g — iA 5 a 6 4 g — iA 3 a63 



(104) 



x g-iA 2 a62g-iA3Q:61 g-iAi7a 6 Og-iA3Q!59g-iAioQ!5Sg-iA3Q!57 

5g-iA 2 a54g-iA3 



V3039 



g -"o — ui-g "'ii^uug -"j-yyg '"iu"jDg «--o"o 

x g-iA 5 a 56 g-iA3Q!55g-iA2a54g-iA3a53g-iA26Q!52g-iA3a5i 
x g— iAi7a 50 g— iA 3 a 4 9 g — iAi O!48 g — iA 3 a 4 7g— iA 5 a 46 g — iA 3 a 4 

">30=43g-iA; 

x e ~ iX " a ™ e - 
x e _i - >V2Q!32 e _i ' >v; 

x g-iA 2 6a26g-«A 



x g — iA 2 Q!44g — iA 3 Q!43g — iA 3 7Q!42g — iA 3 Q!4ig — iA 2 6a40g — iA ; 

»g-i^3a37g-iA 



g-iA 5 0!34giA3Q!33 

g-iA 50 a3 g-iA ; 



■g-iAi a36g-iA; 

>3Q!31g-i^! 



g — 1^13035 g— iA| 

>3Q:29g-iA 



'37C*28g — sA; 



V3027 



^3"25g — iAi7a 2 2g — iA 3 Q!21 g — iAi Q!20g — iA 3 Q!i9 

x g-iA 5 aiSg-iA3Q:i7g-iA65«l6g-iA3 



x g-iA37ai 2 g-iA 3 Q!iig-iA: 

x e - jAl0Q6 e- 



>3Q!l5g— iA 50 Q!i4g— iA 3 ai3 

.26"10g — iA 3 Q!9g — iAi 



g — ./.iragg- iA 3 a 7 
,g-iA3a 5 g-iA 5 Q!4g-iA3a3g-iA2a2g-iA3Q!i 



(0 4 ) 2 sin(# 5 ) 2 sin(0 6 ) 2 sin(0 7 ) 2 sin(0 8 ) 2 

2 * A3 

(1 + 3 cos(2 6» 2 )) sin(6» 3 ) 2 sin^) 2 sin(^5) 

J7I 



lg cos(2 #1) sin(# 2 ) 2 sin(# 3 ) 2 sin 
¥ 2 

2 : " ^ 2 sin(6> 6 ) 2 sin((9 7 ) 2 sin(6> 8 ) 



1 

(1 + 2 cos(2 6> 3 )) sin(6'4) 2 sin(6<5) 2 siring) 2 sin(6< 7 ) 2 sin(0 8 ) 

9 ^/fi 



^A 8 



2V6 

(3 + 5 cos(2 6> 4 )) sin(6» 5 ) 2 sin^e) 2 sin(6» 7 ) 2 sin(6» 8 ) 2 



! A 3 5 



-*Ai 5 



4V10 

(2 + 3 cos(2 6» 5 )) sin(6» 6 ) 2 sin(6» 7 ) 2 sin(6» 8 ) 2 
2VTE 

(5 + 7 cos(2 6> 6 )) sin(<9 7 ) 2 sin(<9 8 ) 2 



l=A24 



4V2T 

(3 + 4 cos(2 6» 7 )) sin(6» 8 ) 2 



=A 48 



(A 



4V7 

3_cos(208)\ 

4 J 



*A63 



*A 80 . 
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From Appendix [b] and equation (^7|) we know that the ranges for the quotient group SU(9)/Zg and the pd parameters 



are 



< a 2i -i < it, < a 2i < |, (1 < i < 36) 

' ^<0 3 <^, (l<j<8), (106) 



,vj+t; " J " 2 

combined with equations ([59]) and (j60|), which yields the Haar measure for SU(9) 

dVsu(9) = sin(2a2) cos(a4) 3 sin(a4) cos(ae) 5 sin(ae) cos(a 8 ) 7 sin(a 8 ) cos(aio) 9 
x sin(aio) cos(ai2) n sin(ai2) cos(ai4) 13 sin(ai4) cos(aig) sin(aig) 
x sin(2ai 8 ) cos(a 2 o) 3 sin(a 2 o) cos(a 2 2) 5 sin(a 2 2) cos(a 2 4) 7 sin(a 2 4) cos(a 2 6) 9 
x sin(a26) cos(a2s) n sin(a2s) cos(a3o) sin(a 30 ) 13 

x sin(2a 3 2) cos(a 34 ) 3 sin(a 34 ) cos(a 36 ) 5 sin(a 36 ) cos(a 38 ) 7 sin(a 38 ) cos(a4o) 9 
x sin(a4o) cos(a42) sin(a42) 1 

x sin(2a44) cos(«46) 3 sin(a46) cos(a4 8 ) 5 sin(a4 8 ) cos(aso) 7 sin(a5o) cos(a52) sin(as2) 9 
x sin(2a54) cos(a 56 ) 3 sin(a 56 ) cos(a 58 ) 5 sin(a 58 ) cos(a 60 ) sin(a 6 o) 7 
x sin(2a 62 ) cos(a64) 3 sin(a 6 4) cos(«66) sin(a66) 5 
x sin(2a6 8 ) cos(a7o) sin(a7o) 3 

x sin(2a72)<ia 8 o . . . doti, (107) 
one can calculate the group volume for the two qutrit system; 



9-1 (9-l)(9 + 2) 

V SU(9) =* 2 7T 2 



fe=l v 7 



(108) 



105345515520000 
which is what one gets with N = 9 from equation (|S4"|). 

XII. CONCLUSIONS/COMMENTS 

The aim of this paper has been to show an explicit Euler angle parametrization for SU(N) as well as the Hilbert 
space for all N by N density matrices. 10 This parametrization also yields a general form for the Haar measure 
for SU(N) as well as confirms Marinov's initial group volume formula for SU(N). It should be noted that such 
a parametrization could be very useful in numerous areas of study, most notably numerical calculations concerning 
entanglement and other quantum information parameters. 11 This parametrization may also allow for an in-depth 
analysis of the convex sets, sub-sets, and overall set boundaries of separable and entangled qubit, qutrit, and N- 
trit systems without having to make any initial restrictions as to the type of parametrization and density matrix in 
question. 
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10 This is but one parametrization that can be realized for a unitary group. For example Dr. Vilenkin and Dr. Klimyk in Representation 
of Lie Groups and Special Functions |19[ have an extensive treatment of the unitary and orthogonal groups using a similar treatment 
as was used here. 

11 This has already occurred for the 5(7(4) case, see for example quant-ph\0203088 and quant -ph\0207181. 
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APPENDIX A: DIAGRAM FORMALISM FOR SU(N) PARAMETRIZATION 

The invocation of "removal of the redundancies" in order to achieve the necessary N 2 — 1 number of parameters 
in our SU(N) parametrization seems to be an arbitrary addition to an otherwise rigorous mathematical discussion. 
This assumed arbitrariness is due, in part, to an underlying construction that is the basis for the formulation of the 
Euler angle parametrization that we have given. The goal of this appendix is to describe that construction and to 
alleviate any remaining sense of arbitrariness in the development of the parametrization. 12 

To begin, let's take a two-component column matrix w given by 



(Al) 



We wish to find the minimal number of operations that rotates, or "mixes" , the two components. Since w is a 2 x 1 
column matrix, our operations must be 2 x 2 square matrices in order to generate a 2 x 1 column matrix as output. 
These 2x2 square matrices are going to be elements of SU(2). Therefore we should look at the well known Euler 
angle parametrization of SU(2) to find these operations. 

Recall that SU(2) is the covering group of 5*0(3), the group of rotations of three dimensional Euclidean space, with 
the following representation: 

Every rotation R € SO (3) can be parameterized by an axis of rotation n and the angle 9 of rotation about 
the axis: R = (n,0). The axis requires two angles, (a,/?), for its specification so three parameters are 
needed to specify a general rotation: 5*0(3) is a three parameter group S. 

Now, any arbitrary rotation R in Euclidean space can be accomplished in three steps, known as Euler rotations, which 
are characterized by three angles, known as Euler angles. When one wants an arbitrary rotation in (x,y,z) space, one 
then does the following: 

1. Rotate through an angle a about the z-axis. 

2. Rotate through an angle (3 about the new y-axis. 13 

3. Rotate through an angle 7 about the new z-axis. 

The rotation matrix R describing these operations is thus given by 

R(a,p,i) = Rz ncw {i)Ry ncw (P)RM), (A2) 

where 

(cos [a] sin [a] 0\ (\ \ / cos [a] sin [a] 0\ 

-sin[a] cos[a] , Ry new (/3) = cos[/3] sin[/3] , R Znevl {-)) = - sin [a] cos[a] . (A3) 
1/ \0 -sin[/3] cos[/3]/ \ 1/ 

Unfortunately at this point we have to make a distinction between body-fixed (the new axes) and space-fixed (the 
original axes) coordinates. For a general rotation, we would like to represent R(a, (3, 7) in terms of space-fixed axes 
^]. ^0| and not body-fixed axes. 14 This is can be achieved in two steps. First one makes the following definitions pi] 

Rih^iP) = RzWRyiftiRMY 1 and i? 2 „„„( 7 ) = Ry n MRz{l){Ry new {0))~ l ■ (A4) 



Then, substituting these definitions back into the body-fixed version of R(a, (3, 7) given in equation flA2j ) we get 

R{a,P,i) = R Zncm {i)R Vncw {(3)R z {a) 

= Rynev {P)RM{Ryn« W OW") 

= Ry„e m (P)Rz(l)Rz(a) 

= R^Ryi^iR^a^R^R^a) 

= R z {a)R v {f3)R z { 1 ) (A5) 



12 Using diagrams to develop a representation of a group is a standard tool in developing group parametrizations. An excellent example 
of this is Dr. Vilcnkin and Dr. Klimyk Representation of Lie Groups and Special Functions section 10.5 [|hJ. 

13 In Classical Mechanics, this rotation is conventionally done about the new x-axis pfj . 

14 This distinction comes from the physical origins of the three-dimensional rotation group: the rotations of a rigid body about a fixed 
point constitute the group which we now know to be SO(3). The terminology is from the physical distinction between the coordinate 
frame "attached" to the rigid body, and the coordinate frame of the "surrounding" system 



25 



where we define the y-axis rotation as 

cos[/3] sin[/ 

Ry(0) = ( 1 | (A6) 
cos[P} / 

and the last step in the derivation was accomplished by noticing that rotations about the same axis commute. 

Thus, one generally does a rotation about the z axis, then about the y axis, then another rotation about the z axis, 
in order to describe the most general rotation of a rigid body. But, if one is rather interested in the "mixing" of the 
components of a Euclidean three- vector describing a point in the rigid body (from the point of view of the space-fixed 
axes), then a rotation about the z axis is equivalent to a "mixing" between the first and second components of the 
three vector, the rotation about the y axis is equivalent to a "mixing" of the first and third components, and the 
second z axis rotation, can be seen to be yet another "mixing" of the first and second components. 15 If we look at the 
second z axis rotation as imparting an overall phase, then the first two rotations become the key "mixing" rotations 
while the third rotation imparts an overall phase to the components of the three- vector. 16 

In conclusion then, we need the relationship between the rotation matrix R(ot,(3, 7) and it's SU(2) counterpart in 
order to accomplish the desired "mixing" of w. This relationship comes from the following mappings between the 
50(3) matrix elements and the SU(2) matrix elements (see for example || |t], |j| for more details) 

'10 



± z:i =* o cosm 



ana 



cosily — sinf^T 



,0 sin [a] cos [a] 
cosf/31 sinf/3 



±1 l ?L by =M 10 

■sm[f] cos[fi] ) \-sin\B] cosf/ 



'cos[7] — sin[7] 0^ 
n[7] cos [7] 
1, 



± ( € q e a ) I sin H cos[ 7 ] I (A7) 



From these mappings we see that for a U S SU(2), an equivalent rotation matrix to R(a, 0, 7) can be written as 17 

e i(9+i>) cog [^ e i(0-VO sin ^] \ 

- e i(-e+^) sin [0] e -i(e+i>) cos [^ )■ \ Aii ) 



If we act upon w with equation (A8) we generate the rotated column matrix W 



_ / a * e^+^ cos[0] + b * e^ e ~^ sin[0] \ 
W _ [b * e -*(»+*) cos [ ] _ a + ^(-0+^) sin ^J ■ l Ay ^ 

Therefore, following the previous discussion concerning SO(3) rotations, a simplified rotated, or "mixed" form of w 
can be seen to be given by 



a * e 10 cos[(/>] + b * e 10 sin[(/>] 
sin[(/>] 



(A10) 



15 We will see a similar discussion when we look at elements of SU{2>) acting on a three component column matrix. 

16 The fact that we are having the second z-axis rotation impart nothing more than an overall phase is due in part to an understanding 
that this geometric approach to the parametrization is not to be taken as as a completely rigorous mathematical derivation, but rather 
as an aid in visualizing the important group elements of the SU (N) Euler angle parametrization (which when combined with the 
correspondingly appropriate exponentiated Cartan subalgebra components at the end of the representation, generates both the correct 
number of parameters for the group, as well as correctly imparts the needed overall phase onto "mixed" column matrix components). 
For example, in the SU(2) parametrization case, the third rotation is considered part of the Cartan subalgebra of the parametrization, 
and can therefore be ignored until later in the paramctrization's development. 

17 The equivalence is that 6=^,0=^, and if) = ^. 
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This reduced SU (2) group operation generating W can be pictorially represented by the following diagram 



(All) 



With this pictorially representation of the group action which "mixes" column matrix components, we can now 
graphically express and decompose iV-component column matrix "mixings." 

We can start using this graphical method by looking at a three component column matrix. Here what we want is 
the minimal number of operations that will effectively "mix" the three component column matrix v 



(A12) 



Comparing with the previous work, we notice that since w is a 3 x 1 column matrix, our operations must be 3 x 3 
square matrices in order to generate the necessary 3x1 column matrix as output. These 3x3 square matrices are 
going to be elements of SU(3). Therefore we should look at the Euler angle parametrization of SU(3), initially given 
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and later in jLlJ, to find these operations. 
Before we address the full SU(3) group, the necessary rotations of v can be expressed by extending our diagram 
representation of the SU(2) rotations in the following way 



(A13) 



a 




a 
b 




a 
h 




a 




b 




then ( 


and 
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c 












c 





where 



and 



a 




a 




a 
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)- 




then 


c 
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(A14) 



Since we wish to rotate, or "mix" , components a and b and components b and c, we can see that this can be done in 
stages; first by mixing components a and b then by mixing a and b and a and c together. Therefore, what is needed 
to accomplish this scheme is the group element of SU(3) that will rotate components a and c, without rotating b 
(although an overall phase on b is acceptable). This action is represented by the following SU(3) element 



cos[£] sin[£] 
= I 1 
— sin[C] cos[C] 



(A15) 



With equation ( A15) and following the methodology from the SU(2) work, we can associate the previously diagrammed 
SU(3) rotation element with the following group element from SU(3) 



(A16) 



Thus the SU(3) rotations, diagrammed in equation ( A13 ), are seen to be equivalent to the following group action 19 



a 








a 










a 








a 


b 




b 








b 


c 








c 








= e 




i\2 


<t>2 e 


jA; 



(A17) 



The reason we have a A3 and not a As in the first group element is beyond the scope of this Appendix, but is addressed in the following 
Appendices. 

In this representation we are using the general extension of the fact that the "product" of two rotations Ri and R2, denoted R1R2 is 
the transformation resulting from acting first with R2, then with R\ [Jl3|. 
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The extension of this work to "mixing" four component, and thus, n-component column matrices is straightforward. 
For a four component column matrix t 



b 



(A18) 



we schematically represent the necessary rotations as 



a 
b 




a 
b 
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then ^ 
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(A19) 



where 
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then 




then 
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(A20) 



Following the pattern of the previous work, the rotation between components a and d is achieved by acting upon t 
with a 4 x 4 matrix which is going to be a SU (4) group element 



( cos[x] sin[x] s 

10 

10 

V-sin[x] cos[x]y 



(A21) 



Mathematically then, the above "mixing" scheme can be seen to be given by the following group action 
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'a 




a 




a 




a 






b 




c 




A 




b 




c 




b 
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i\2 


03 g 









a 
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A 

(A22) 

In general then, our graphical representation of "mixing" column matrix components comes from both Biedenharn's 
work on angular momentum |l3[ and from a conceptual extension of the Clebsch-Gordon coefficient construction for 
spin i particles Recall that for spin i particles, the only part of the SU(2) rotation operator which mixes different 
quantum number m values in the Clebsch-Gordon coefficients is the middle rotation; the hrst and third add nothing 
but an overall phase. As we have seen, that middle rotation corresponds to the <T2 Pauli Matrix for SU(2) and to the 
\(N-i) 2 +i matrix for SU(N), N > 2. It is these Lie algebra components that when exponentiated, "mix" the column 
matrix components in question. One then only needs to "attach" (matrix multiply) the appropriate exponentiation 
of the Cartan subalgebra ( 0-3 fo r SU(2), A3 and As for SU ( 3), and A3, As and A15 for SU (4 )) at then end of the group 
action (given in equation (All) for SU(2), equation (A17) for SU(3), and equation (A22) for S'C/(4)) to achieve the 
necessary overall phase for the column matrix components. The generalization of this methodology to higher SU(N) 
groups should now be quite apparent. 



See previous footnote. 
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APPENDIX B: INVARIANT VOLUME ELEMENT NORMALIZATION CALCULATIONS 

Before integrating dVsu(N) we need some group theory. We begin with a digression concerning the center of a 
group jj], |22j . If S is a subset of a group G, then the centralizer, Cq(S) of S in G is defined by 

C(S) = C G (S) = {xeG\ if s G S then xs = sx}. (Bl) 

For example, if S = {y}, C(y) will be used instead of C({y}). Next, the centralizer of G in G is called the center of 
G and is denoted by Z{G) or Z . 

Z(G) = Z = \z G G I zx = xz for all x G G} 

= W <B2) 

Another way of writing this is 

Z{G) = n{C(x) | x G G} 

= {z | if x E G then z G C(ar)}. ( ' B3 ^ 

In other words, the center is the set of all elements z that commutes with all other elements in the group. Finally, 
the commutator [x, y] of two elements x and y of a group G is given by the equation 

[x,y] = x~ 1 y~ 1 xy. (B4) 

Now what we want to find is the number of elements in the center of SU(N) for N — 2, 3,4 and so on. Begin by 
defining the following 

Z n = cyclic group of order n = Z„ Z(SU(N)). (B5) 

Therefore, the set of all matrices which comprise the center of SU(N), Z (SU (N)) , is congruent to Z^ since wc know 
that if G is a finite linear group over a field F, then the set of matrices of the form Ec g <?, where g G G and c 9 G P 1 , 
forms an algebra (in fact, a ring) || ^2|. For example, for SU(2) we would have 

Z 2 ={x e S?7(2) | G Zi for all y G 5?7(2)}, 

[a?,y] =wl3, (B6) 
Z x ={1 2 }. 

This would be the set of all 2 by 2 matrix elements such that the commutator relationship would yield the identity 
matrix multiplied by some non-zero coefficient. In general this can be written as 

Z N ={x G SU(N) | [as, y] G Z x for all y G SU(N)}, 

zi={i»}. (B7) 

This is similar to the result from 0, that shows that the center of the general linear group of real matrices, GLpf(JHt), 
is the group of scalar matrices, that is, those of the from ujI, where I is the identity element of the group and u> is 
some multiplicative constant. For SU(N), tut is an N th root of unity. 

To begin our actual search for the normalization constant for our invariant group element, we first again look at 
the group SU{2). For this group, every element can be written as 

% i) < B8 > 

where \a\ 2 + \b\ 2 = 1. Again, following [j7| we can make the following parametrization 

a = Vi- m, 

b = y 3 ~ m, (B9) 

i = vl + yl + vi + vl 

The elements (1,0,0,0) and (—1,0,0,0) are anti-podal points, or polar points if one pictures the group as a three- 
dimensional unit sphere in a 4-dimensional space parameterized by y, and thus comprise the elements for the center 
group of SU(2) (i.e. ±12). Therefore the center of SU(2), Z2, is comprised of two elements; ±12- 
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Now, in our parametrization, the general SU (2) elements are given by 

D(u,v,£) = e^e^V^, , x 

' (BIO) 
dVsu(2) = sm(2v) dfidvd£, 

with corresponding ranges 

0</i,£<7T, (Bll) 

0<i/<|. (B12) 

Integrating over the volume element cZVgfy^) with the above ranges yields the volume of the group SU(2)/Z2- In 
other words, the SU(2) group modulo its center Z-i. In order to get the full volume of the SU(2) group, all ones need 
to do is multiply the volume of SU(2)/Z2 by the number of identified center elements; in this case two. 

This process can be extended to the SU(3) and SU(4) parametrizations. For SU(3) ]9], [l(], |l^, [l7|] (here recast as 
a component of the SU (4) parametrization derived in p| ) 

SU(3) = e tX3a7 e tX2a *e tX3a9 e lX5aw D( ail , a 12 , a 13 )e lXsai \ (B13) 

Now, we get an initial factor of two from the D(an, a\2, dxs) component. We shall now proves that we get another 
factor of two from the e lA 3 Q 9 e ' A5Ql0 component as well. 

From the commutation relations of the elements of the Lie algebra of SU{3) (see || for details) we see that 
{A 3 , A 4 , A 5 , A 8 } form a closed subalgebra SU(2) x U(l). 21 

[A 3 , A 4 ] =1X5, 
[A3, A 5 ] = - iA 4 , 
[A3, A 8 ] =0, 

[A 4) A 5 ] =i(A 3 + V3A 8 ), 
[A 4 , Ag] = - «v3A 5 , 

[A 5 ,A 8 ] =iV3A 4 . (B14) 

Observation of the four A matrices with respect to the Pauli spin matrices of SU(2) shows that A 4 is the SU(3) 
analogue of <7i, A5 is the SU(3) analogue of 02 and As is the SU(3) analogue of 173 



CTl 



1 

1 




1 
1 

-2, 



Thus one may use either {A3, A5} or {A3, A5, As} to generate an SU(2) subgroup of SU(3). The volume of this SU(2) 
subgroup of SU(3) must be equal to the volume of the general SU{2) group; 2vt 2 . If we demand that any element of 
the SU(2) subgroup of SU(3) have similar ranges as its SU{2) analogue 22 , then a multiplicative factor of 2 is required 
for the e iX3<xs e i\sa w com p 0n ent. 23 



21 Georgi [ |l 2| has stated that A2, A5 , and A7 generate an £C/(2) subalgebra of SU (3). This fact can be seen in the commutator relationships 
between these three A matrices contained in 1l| or in M. 

22 This requires a normalization factor of -j= on the maximal range of Ag that is explained by the removal of the Z3 elements of SU (3) . 

23 When calculating this volume element, it is important to remember that the closed subalgebra being used is SU(2) X U(l) and therefore 
the integrated kernel, be it derived either from e 1 ^ 3 " e^s^e 1 ^ 37 or e tX3a e tX5 ^e lX8 " ! , will require contributions from both the SU(2) and 
U(l) elements. 
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Finally, SU(3) has a Z 3 whose elements have the generic form 

| . (B16) 

\o o 7? 1 - 1 %7 

where 

»7i = % 3 = 1. (B17) 
Solving for rji and 772 yields the following elements for Z 3 
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(B18) 

which are the three cube roots of unity. Combining these SU (3) center elements, a total of three, with the 2 factors 
of 2 from the previous discussion, yields a total multiplication factor of 12. The volume of SU(3) is then 

V su{3) = 2*2*3*V(SU(3)/Z 3 ) 

= V3ir 5 (B19) 

using the ranges given above for the general SU(2) elements, combined with < au < -^=. Explicitly: 

< a 7 , ctg,aii,cti3 < 7r, 
< as, aio, ai2 < — , 

< a X4 < 4=. (B20) 
v3 

These are modifications of || |l0|, |l(| [l7], |23| and take into account the updated Marinov group volume values j|. 

For SU(4) the process is similar to that used for SU(3), but now with two SU(2) subgroups to worry about. For 
a U G SU(4), the derivation of which can be found in El, we see that, 

U = e lA3Ql e iA2Q2 e iA3Q3 e iA5Q4 e lA3Q5 e lAl0Q6 [S';7(3)]e lAl5Q!15 . (B21) 

Here, the two SU (2) subalgebras in SU(4) that we are concerned with are {A3, A4, A5, As, A15} and {A3, Ag, Aio, As, A15}. 
Both of these SU(2) X 17(1) x U(l) subalgebras are represented in the parametrization of SU(4) as SU(2) subgroup 
elements, e iA3Q3 e iA5Q4 and e iA3Q 5 e 'A 10 a 6 _ We can see that ^ Q is the 5^7(4) ana logue of <j 2 24 and A i5 is the SU(A) 
analogue to <r 3 25 . The demand that all SU{2) subgroups of SU(A) must have a volume equal to 27r 2 is equivalent to 
having the parameters of the associated elements of the SU (2) subgroup run through similar ranges as their SU (2) 
analogues. 26 As with SU(3), this restriction yields an overall multiplicative factor of 4 from these two elements. 27 
Recalling that the SU(3) element yields a multiplicative factor of 12, all that remains is to determine the multiplicative 
factor equivalent to the identification of the SU(4) center, Z4. 

The elements of the center of SU(4) are similar in form to the ones from SU(3); 

(B22) 
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We have already discussed A5 in the previous section on SU(3). 
It is the SU (4) Cartan subalgebra element. 

This requires a normalization factor of -4= on the maximal range of A15 that is explained by the removal of the Z4 elements of 5(7(4). 
When calculating these volume elements, it is important to remember that the closed subalgebra being used is SU(2) X U(l) X U(l) 
and therefore, as in the SU(3) case, the integrated kernels will require contributions from appropriate Cartan subalgebra elements. 
For example, the e*^ 3 " 3 e 1 ^ 5 ™ 4 component is an SU(2) sub-element of the parametrization of 5C/(4), but in creating its corresponding 
SU(2) subgroup volume kernel (see the SU(3) discussion), one must remember that it is a SU(2) C SU(3) C SU(A) and therefore the 
kernel only requires contributions from the A3 and \g components. On the other hand, the e lA3Q5 e lAl0 " 6 element corresponds to a 
SU(2) C SU(A) and therefore, the volume kernel will require contributions from all three Cartan subalgebra elements of SU(A). 
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where 



vi 



4 



4 



i. 



(B23) 



Solving yields the 4 roots of unity: ±14 and ± H4, where I4 is the 4x4 identity matrix. So we can see that Z4 gives 
another factor of 4, which, when combined with the factor of 4 from the two SU(2) subgroups, and the factor of 12 
from the SU(3) elements, gives a total multiplicative factor of 192. 

Thus, when we integrate the ST/ (4) differential volume element with the ranges given previously for the general 
SU(2) and S£/(3) elements combined with the appropriate range for the A15 component (all combined below) 

< ai,a 3 ,a 5 ,a 7 , ag, a n , «13 < tt, 
< a%, 0-4, a 6 , a s , a w , a i2 < — , 



0<au<^, 



< ai 5 < 



y/6' 



(B24) 



we get 



Vf 



SU(i) 



2*2*2*2*3*4* V(SU(4)/Z 4 ) 
V2n 9 



(B25) 



This calculated volume for S/7(4), the derivation of which can be found in |Q agrees with that from Marinov j3|. 

From this work, it is plain to see that in general, the ranges for the A3, A2 analogues (recall equation (p0|)), and the 
remaining Cartan subalgebra components of the parametrization will take the following general form 

< a(A 3 ) < 7T, 



< a(A( fc _ 1 )2 + i) < 



< a( Cartan subalgebra elements ) < tt. 



k(k-l) J 



(B26) 



for 2 < k < N. Also, it is apparent that the elements of the center of SU (N) will have the form 
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V 0' 
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where 



N jv jv \ 1 

Vi =V2 =V 3 =■■■ = Vn-i = 1) 



Vi 1 il2 S3 1 ■ 



N 



\ 



(B27) 



(B28) 



thus yielding N elements in Z{SU(N)) and a corresponding multiplicative constant of N in the calculation of the 
group volume. Observation of the previous calculations of the invariant volume element for SC/(2), SJ7(3), and 
SU(A) also indicates that the Euler angle parametrization for SU(N) yields (N — 2) SU(2) subgroups that require 
multiplication by two in order to satisfy the 2ir 2 general SU(2) volume requirement. Therefore, if one defines the 
ranges given in equation (B2£) as V , we see that the invariant volume element for SU(N) can be written as 



SU{N) 



= fijv * 



dV, 



SU(N) 



(B29) 



An observant reader will notice that we have counted A3 twice in the above ranges. The fact of the matter is that we wanted to distinguish 
between the first and second A3 elements in the fundamental SU(2) parametrization. The first is given as the o(A3) component while 
the second is regarded as the SU(2) Cartan subalgebra component. The reason for this distinction will be made clearer in Appendix O. 
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where 

■>iV-2 



n N = 2"- z *N*n N ^ 1 (B30) 



since the differential volume element dVsu(N), given in (|37]), shows a reliance on the differential volume element of 
SU(N — 1) and therefore on fljy-i. 
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APPENDIX C: MODIFIED PARAMETER RANGES FOR GROUP COVERING 

In order to be complete, we list the modifications to the ranges given in Appendix B] that affect a covering of SU(2), 
SU(3), SU(4), and SU(N) in general without jeopardizing the calculated group volumes. 
To begin, in our parametrization, the general SU(2) elements are given by 

D(jjL,u,Z) = e »»e iX > I 'e iX *t, 
dVsu(2) — sia(2v) dfidud^, 
with the corresponding ranges for the volume of SU(2)/Z2 given as 

0</U,£<7T, 

0<v<~. (C2) 

In order to generate a covering of SU(2), the £ parameter must be modified to take into account the uniqueness of 
the two central group elements, ±l2 5 under spinor transformations. This modification is straightforward enough; f's 
range is multiplied by the number of central group elements in SU(2). The new ranges are thus 

< ft < it, 

7T 

0<u<-, 

< £ < 2tt. (C3) 

These ranges yield both a covering of SU(2), as well as the correct group volume for SU(2). 29 

One can also see this by looking at the values of the finite group elements under both sets of ranges. To do this, 
we first partition D(/x, v, £) as 

D({p, v), = ( e « A ^e lA2 ")e lA3e , (C4) 

where 



J\ 2 V = 



e a 3 £ = 



e l » 
e"^ 

cos{v) sin{v)\ - C5 . 
—sin(v) cos{v) J ' ^ ' 



e" ,; « 



Then looking at < £ < 2ir first we see the following pattern emerge 

e— — r-o 1 _°o— 



C-2tt=^(q ^)=^1 2 ^=C = 0. Repeat ( C'(i ) 



Next, we look at < \i < n and < v < ? 



M = , i/ = =>- 1 2 , 

fj, = TT , V = =>■ -1 2 , 

M = 7 r ! ,= |^(?- 1 V (C7) 



One may interchange and £'s ranges without altering either the volume calculation, or the final orientation of a two-vector under 
operation by D. This interchange is beneficial when looking at Euler paramctrizations beyond SU(2). 
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Combining these two results yields the following 12 matrices for D((/j, v), £) 

D((0,0),0)=( I ?V D((7r,0),0) ' ' " 



C((0,0),7r) = -(J J), D((7r,0),7r)=(j J 

D((0,0),27r)=(j J), D((7r,0),27r)=-(j J 

£>((0, f),0) = (_°x J). ^(t>|) ) 0)=(J 

D((0,f),7r) = (5 "o 1 ), D(( 7 r,|), 7 r)=(_° 1 J 

D((0, f ), 2vr) = ( A ) ' ^ I } ' 2?r) = ( 1 ~0 ) ' (C8) 
We can see that in the above 12 matrices we have repeated four fundamental forms, each three times 

-la =*> 3, 3, (C9) 

when we have run £ from to 2n. But, if we remove the repeating I2 term in the original £ calculations above - in 
essence forcing £ to range from to tt instead of 2n - we see that we then get the same four fundamental forms, but 
now each only repeated twice 

-h => 2, (I - 1 ) 2. (CIO) 

The greatest common divisor of the above list is, obviously, 2, which not only corresponds to the amount that the 
range of the £ parameter was divided by, but also to the multiplicative factor of 2 that is required in the calculation 
of the invariant volume element when using the quotient group ranges. This may seem trivial, but let us now look at 
SU(3). 

For 51/(3), here given as a component of the SU(4) parametrization, we know we have one SU(2) subgroup 
component (from Appendix |b|) as well as the SU(2) contribution _D(/x, z/, £) here rewritten in terms of the SU(3) 
parameters an, ai2, and 0*13: 

517(3) = e lA3a V A2Q V A3Q V A5Ql0 i}(aii, a 12 , a 13 )e lXsai4 . (Cll) 

Therefore the ranges of ag and aia should be modified just as £'s was done in the previous discussion for SU(2). 30 
Remembering the discussion in Appendix |B| concerning the central group of SU(S), we can deduce that cum's ranges 
should be multiplied by a factor of 3. This yields the following, corrected, ranges for SU(3) 31 

< «7, an < 7r, 

7T 

< a s , atio, 012 < — , 



< a 9 ,a 13 < 2ir, 

< a u < VSn. (C12) 



Recall the previous footnote on the interchange of the first and third rornpo r ien t ranges in a SU (2) parametrization. 

Earlier representations of these ranges for SU(3), for example in H, llOl llsL h/H, E3] , were incorrect in that they failed to take into account 

the updated SU (N) volume formula in B . 
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These ranges yield both a covering of SU(3), as well as the correct group volume for SU(3). 
For aUe 5£/(4), we have two SU(2) subgroup components 

U _ e i\3<X-L e %\20C2 e %\30C3 gi^SOCl e i\30ts e i\-L0Ot6 [gf/^)] e*^ 15 " 15 (C13) 

As with the SU{2) subgroup ranges in 5?7(3), the ranges for a$ and a§ each get multiplied by 2 and ais's ranges get 
multiplied by 4 (the number of SU (4) central group elements) . The remaining ranges are either held the same, or 
modified in the case of the 517(3) element; 

< ai, 07, an < tt, 

7T 

< a 2 , 0!4, "6, "8, «10, "12 < 2", 
< a 3 , a-5, a g , ai3 < 2-7T, 
< ai4 < \/37r, 

< aie < 2^77. (C14) 

These ranges yield both a covering of SU(4), as well as the correct group volume for SU(4). 

In general we can see that by looking at SU(N)/Zm not only can we arrive at a parametrization of SU(N) with a 
logically derivable set of ranges that gives the correct group volume, but we can also show how those ranges can be 
modified to cover the entire group as well without any arbitrariness in assigning values to the parameters. 
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